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We investigate the electronic transport through a suspended carbon-nanotube quantum dot. In 
the presence of a magnetic field perpendicular to the nanotube and a nearby metallic gate, two forces 
act on the electrons: the Laplace and the electrostatic force. They both induce coupling between 
the electrons and the mechanical transverse oscillation modes. We find that the difference between 
the two mechanisms appears in the cotunneling current. 

PACS numbers: 71.38.-k, 73.63.-b, 73.63.Kv, 85.85.+j. 



I. INTRODUCTION 

Displacements of conducting nano-mechanical systems 
are conveniently registered by the electron transport 
measurements. The electron transport is sensitive to 
the mechanical motion mostly due to the electrostatic 
interaction, which allowed one to detect, e.g. the mo- 
tion of tiny single systems, like Ceo molecules^ or car- 
bon nanotubesi^ The latter are particularly promising 
systems for many applications. High mechanical quality 
factors of the order of 10^ have been reported for sus- 
pended carbon-nanotube oscillatorsi^ At the same time, 
strong mechanical coupling to electronic transport, both 
in the bendingii^ and breathing modesj^ii have also been 
demonstrated. The effect of the electrostatic coupling 
to the mechanical oscillations (also known as polaronic 
coupling) on transport has been investigated by several 
authors from the theoretical point of view<^"— A suppres- 
sion of the current at low bias voltage, known as Frank- 
Condon blockade, has been predicte d^^i^^ and observed^ 
in the quantum regime, for high frequency molecular os- 
cillators. Similar effects in the classical limit (low fre- 
quency oscillators) have also been discussedJi"— 

Much less investigated is the effect of magnetic field, 
see Fig. [1] on the coupling between electronic transport 
and mechanical oscillations (called magneto- elastic cou- 
pling hereafter) . This coupling is at the origin of the mag- 
netomotive effect, which has been used to activate and 
detect the classical motion of micro and nano-mechanical 
resonators since a long time. However its manifestation 
in the quantum regime has been addressed only recently 
for a suspended carbon nanotube forming a quantum 
doti^"— A particularly appealing interpretation of the 
magnetoresistance in this system has been put forward 
in Ref. [l^: In the tunneling limit the electrons travers- 
ing the nanotube interfere over the different paths gener- 
ated by its zero point oscillations. The reduction of the 
current can then be seen as a joint manifestation of the 
Aharonov-Bohm effect and the quantum fluctuations of 




FIG. 1: Schematic picture of the system: a voltage biased 
suspended carbon nanotube in presence of a transverse mag- 
netic field B and of a transverse electric field E created by 
applying a gate voltage. 



the mechanical degree of freedom. In the resonant case 
the magnetoresistance has been calculated as an expan- 
sion up to the second order in the small electromechanical 
coupling. It has been found that multiple interferences 
increase the size of the effect parametrically<2^ 

In any realistic device both the electrostatic and mag- 
netic effects are present. The aim of the present paper is 
to discuss the interplay of these two interactions. Specif- 
ically the system considered is shown schematically in 
Fig. [T] It consists in a suspended carbon nanotube form- 
ing a quantum dot contacted electrically to two exter- 
nal leads. We focus on the coupling of the fundamental 
bending mode of the nanotube to the electronic transport 
in presence of both a strong magnetic field perpendicu- 
lar to the nanotube and a metallic gate. They induce 
a magneto-elastic and a polaronic coupling, respectively. 
We consider a high temperature regime when electron 
tunneling events between the leads can be described in 
terms of rate equations i2L2^ On the other hand, the me- 
chanical mode is in the low-temperature limit and, for 
simplicity, we will further assume that its coupling to the 
environment is sufficiently strong to keep the oscillator in 
equilibrium. 

The theory that is formulated in the reminder of the 
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paper allows to explore the purely magnetic problem 
outside the limits of previous works)^"— by investigat- 
ing transport at bias voltages larger than the oscillator's 
energy quantum. The same approach is then used to 
study the interplay between the magneto-elastic and po- 
laronic couplings. In the sequential tunneling regime, we 
find that the current depends on a single coupling con- 
stant which is the root mean square of the polaronic and 
magneto-elastic coupling constants, respectively. The 
presence of a magnetic field at this order thus only renor- 
malizes the effective interaction and introduces a Zeeman 
splitting of the electronic levels. 

The difference between the two interactions appears in 
the cotunneling regime (when electrons tunnel directly 
between the leads through virtual states in the dot). 
When the single energy level in the dot is close or in the 
conducting window between the left and right chemical 
potential of the leads (resonant transport) cotunneling 
events induce a small correction to the current. How- 
ever, they can dominate transport away from resonance. 
It is thus in that regime that the difference between the 
polaronic and magneto-elastic couplings becomes clear- 
cut and could be detected experimentally. 

The paper is organized as follows. Section |TT] is dedi- 
cated to the microscopic derivation of the model Hamil- 
tonian used in the reminder of the paper. Section IIIII 
introduces the formalism used to tackle the problem. In 
Section IIVI we discuss the behavior of the current as ob- 
tained in the sequential tunneling approximation and the 
consequences of the Zeeman splitting. Section |V] consid- 
ers the contribution of the cotunneling events, discussing 
in particular the difference between the polaronic and 
magneto-elastic couplings. Section IVII gives our conclu- 
sions. Part of the material is presented in the appendices. 



II. MODEL HAMILTONIAN 

We derive the Haniiltonian for a suspended nanowire 
weakly electronically coupled to leads in presence of po- 
laronic and magneto-elastic effects in Sees. Ill AIIII Cl see 
Eqs. (pn | - (PT|) . We estimate the amplitudes of these ef- 
fects for a quantum dot formed with a suspended carbon 
nanotube in Sec. Ill Dl 



A. Magneto-elastic coupling 

In this subsection we present a detailed microscopic 
derivation of the magneto-clastic interaction. It extends 
the model proposed in the previous worksi^^^ 



1. Electron in a vibrating nanowire 

The Hamiltonian that describes an electron propagat- 
ing in a ballistic, suspended nanowire in presence of a 



transverse magnetic field (units with h = ki 
used throughout the paper) is 
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^ = ^ (p + eA{r)f + Vix, y - u{x),z) + F.tdng. (1) 

Here, p = {pxTPy,Pz) and r = (x,y,z) are the momen- 
tum and position of the electron of mass m and charge 
— e, A = {—By ,0,0) is the vector potential associated 
with the magnetic field B = (0, 0, B). The potential 

V{x, y - u{x), z) = Vb{x) + Vconfiy - u{x),z) (2) 

contains two terms. The first one describes the tunnel 
barriers between the nanowire and the leads, 



V^,{x) = Vi5{x) + Vr5{x - L), 



(3) 



where L is the nanowire's length. The second one is 
the confining potential Vconf inside the nanowire, it de- 
pends on the deflection u{x) of the nanowire along the 
y-direction. 

The elastic bending of the nanowire in the (x, ?/)-plane 
is described within the elastic string theory^S with the 
term 



string 



ds 



TT 



EI (d'^u 



(4) 



where p is the mass per unit length of the string, E, its 
elastic Young's modulus, and /, its area momentum of 
inertia, the deflection field and momentum density oper- 
ators of the string are conjugated [u{s), 7r(s')] = i5{s—s'). 
For a doubly clamped nanowire, the boundary conditions 
are u{s) = dsu{s) = at s = 0, L. 

The classical mechanical eigenmodes of the nanowire 
arc characterized by a frequency w„ - the lowest mode 
has frequency loq sa 22A{EI / pL^y^^ - and displace- 
ment field /n(s) (integer n > 0) with normalization 

(1/L) dsfn{s)^ = I. Decomposing the deflection field 
operator onto these modes. 



, Xnfn{s), 



the Hamiltonian ^ then becomes: 



E 



/ p2 

\2M 



MujiXi 



(5) 



(6) 



where Xn and P„ are the effective positions and momenta 
of the modes, with conjugation relation [Xn, Pm] = iSnm, 
and M = pL is the mass of the nanowire. 

We can reasonably expect that the confining potential 
enforces the replacement y — > u{x) in Eq. ([T}, assum- 
ing that the nanowire stands along a line with y = in 
classical equilibrium. The one-dimensional propagation 
of the electron and its coupling to the vibrations would 
then be described by 



n 



2m 



[p,~eBuix)]^ + Vbix)+H, 



string ■ 



(7) 
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A unitary transformation H 



ne 



-iS 



where S 



-eBU{x), with d^U{x) = u{x) would then yield: 



2M 



(8) 



where Ho = pl/{2m) + Vh{x) and F„(x) = dsfn{s), 
which will be the starting point of the next section. 

Appendix |X] is devoted to a rigorous derivation of 
Eq. (IS]) starting from Eq. ([T]) and taking into account 
the transverse spreading of the electron's wavefunction 
in the nanowire. 



2. Single level regime 

Here we simplify further the Hamiltonian ([8]) in 
the regime of a single resonant electronic level in the 
nanowire. 

Indeed, the term Hq in Eq. ([8]) describes an electron 
propagating through a double-barrier system. For strong 
tunnel barriers, quasi-localized states are formed in the 
central region. Due to the coupling with the leads in the 
outer regions, they decay with a lifetime l/F. Assuming 
a large energy level spacing A 3> F between the quasi- 
localized states, we can focus on a single level with energy 
Ed in the central part of the nanowire, provided that the 
chemical potential jj, in the leads is also close to the single 
level energy, \fi - ed\ < A. 

In this regime, the single-particle Hamiltonian en- 
tering Eq. ([5]) can be written equivalently 



Hi 



where 



Ht 



^(t.4,rf + h.c.) 



(9) 



(10) 

(11) 
(12) 



describe the leads, the dot, and their tunnel coupling, 
respectively. Here, Ckv {i' — l,r) are the annihilation 
operators for the electron propagating states in the left 
(I) and right (r) leads, labeled with momenta k, and 
energies = k'^/2m; d is the annihilation operator in 
a single level state of the nanotube, with energy Sd = 
7r^p^/(2mL^) {p is a strictly positive integer). The wave- 
function associated with this level is \/2/L sin(7rpx/i). 
The tunnel coupling between this state and the contin- 
uum of states in the leads is accounted with the tun- 
nel matrix elements t^. Each lead contributes to the 
single- level energy broadening = 2t: N i,\t , respec- 
tively, where Nu are the densities of states (per spin) at 
the Fermi level in the leads (F = F; + F^ is the total 
braodening). They can be related with the parameters 
in Eq. ^ through the relations F^, = 2-^-^ /\{^/V^kpY 



aX V = l,r, where A w n'^p/mL'^ (at large p) and kp is 
the Fermi wavevector. 

We use the basis of electron states defined above to 
write the resolution of the identity (for the single-particle 
problem considered here) as 1 = ni + Ur + Ud, where 
Ui, = c^^Ckiy and Ud ~ d^'d count whether the electron 
is in the left lead, right lead, or the dot, respectively. 
Inserting it on the right and left sides of Eq. - and 
also assuming that A 3> for the relevant oscillator's 
eigenmodes, we find that the Hamiltonian coupling the 
states with energies close to the Fermi level reads 

n = no + ^{MulXl/2 + K,,ini + K^rTir + K^dnd), 

n 

(13) 

with 



(14a) 
(14b) 





[Pn 
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K — 


[Pn 


^eBF^{L)f 




2M 


Knd = 


[Pn 






2M 



2M 



(14c) 



and 



{F:) = \j\xsin^[^)F^{x). (15) 
for a = 1,2. 

In order to eliminate the i3-dependcnt terms in Eqs. 
(|14p . we perform another unitary transformation % — > 
e'^'ne-^^', with 

S' = -eBY,[riiFn{0) + n,F„(L) + {FM X„, (16) 

n 

and we find 



n 



HL+nd + Hr- 



Here, 



where 



iid = edd^d, 



ed = ed + ^{{F;^)-{F.ny') 



(17) 
(18) 

(19) 



includes a diamagnetic shift for the single level, and 
= oxp{ieB ^ X„[(F„) - F{0)])cld + h.c.} 

k ri 

+ Y,{tr cxp(ieB XnKFn) - F{L)])cld + h.c.}. 

k n 

(20) 

Introducing the bosonic annihilation operators 5„ for 
the oscillator's modes (phonons), with Xn = Xno{bn+blJ 
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{Xno = l/\/2MuJn is the amplitude of zero-point motion 
in the eigenmode n), we write 



-^string — ^nbnb„ 



(21) 



We also notice that the eigcnmodes of the doubly 
clamped nanotube have a definite parity and, thus, 
Fn{L) - (Fn) = {-y'[{F.n) - F„(0)]. We then rewrite 
Eq. ^ as 



= 5Z{*-cxp[^^(&„ + 6t)0,„]cLrf + h.c.}(22) 



where </>/„ = = and 



(23) 



is the effective magnetic flux, in units of the flux quan- 
tum $0 — 27r/e, that characterizes the coupling between 
the electron and the eigenmode n of the oscillator. Inte- 
gration by parts allows simplifying Eq. (j23p into: 



ttBX„oL dx 



(24) 



at even n, while a similar formula with an extra factor 
1 — 2a;/i -|- sin(27rp2;/L)/(7rp) in the integrand holds at 
odd n. 

Equation ([T7)) together with is the final result of 
this section. It generalizes the Hamiltonian derived in 
Ref. H^that only considered the lowest eigenmode. [Note 
that different conditions are obtained to justify its valid- 
ity (see Appendix|A| and that the phase (|24|) is smaller by 
a factor 2tt here.] The phase was interpreted there as 
an Aharonov-Bohm phase picked by the electrons when 
they cross the tunnel barriers, and which depends on the 
deflection of the nanowire. 



B. Polaron coupling 

In this subsection we generalize the model to include 
the electrostatic coupling and the presence of the spin 
degree of freedom. We begin by adding an electrostatic 
contribution to the Hamiltonian ((22|>^ 



Hr 



eVnN, 



(25) 



where Cg is the gate capacitance (we assume that the 
lead capacitances are negligible with respect to the gate 
capacitance), Vg is the gate voltage, and N is the opera- 
tor for the total number of electrons on the dot. Since we 
assume that only a single level contributes to transport, 
it is possible to write N = Nq + rid, where Nq « CgVg/e 
is the number of the excess electrons in the filled levels, 
and Ud is the number of electrons in the relevant level. 



For small displacement of the oscillator we can express 
Cg[u[s)\ in terms of A"„: 



Cg[u{s)] = C,[0] 1 



n J 



where 



ds 1 dCg 
L Cg[Q\ du 



(26) 



(27) 



Substituting the expansion into (P5|) one obtains 
(apart from a constant term): 



Hr 



a„X„ -I- U^nd-\nd]. + {5s - eVg) na + Hpoi , 

n 

(28) 

where Uda = d^d^, Ud = ridt + "^di, Uoo = e^/Cg[Q\, 
5e = Uoo{2Na + l)/2, and 



Hpo\ = ^ A„w„(fe„ + bl' 



nd , 



(29) 



with A„ — a„SeXno/uJn- The first term of Eq. (|28p 
gives a shift in the equilibrium position of the oscilla- 
tor and can be discarded. The second term gives the 
intradot Coulomb repulsion. We will assume that Uac is 
the largest energy scale of the problem, allowing to ne- 
glect the contribution of the double occupancy state [the 
term proportional to Uoo times X„ can be neglected for 
the same reason and it is not shown in Eq. ([28))]. The 
third term can be included in the definition of Sd since 
it gives only a renormalization of the energy of the level. 
The last term, defined in Eq. (|29p , is the seeked polaronic 
coupling, which models the interaction between the os- 
cillation modes and the charge fluctuations in the dot. 

The expression for the polaronic coupling constant 
can be calculated by using a model of a distributed 
capacitancei^i^ 



CgHs)] 



L 

ds 



27reo 



(30) 



/o arccosh([(/i — w(s)]/_R^) ' 

where h is the distance of the nanotube from the sub- 
strate, R± its radius, and eq is the vacuum permittivity. 
Substituting Eq. pop in the definition of A„ and calcu- 
lating the integral for u = and R± ^ h one obtains: 



A. 



eCgmVg\ X^o [''ds 
2'ireohL ut^ 



-Us) 



(31) 



for Cg[0]Vg > e. 

Including the electrostatic coupling and the spin de- 
grees of freedom to the Hamiltonian (p2|) we obtain: 



H — iJstring + Hl + Hd 



Hj 



(32) 



where the term -ffstring remains unchanged, see Eq. (|21l) . 
The inclusion of the spin and the chemical potentials 
gives instead 



(33) 
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where c^ka is the annihilation operator for a state with 
momentum k and spin a =t,i in lead v = l,r, = 
£k — fJ-ui and the difference of the chemical potentials 
l-ti — IJt = e.V is related to the bias voltage V of the 
junction. The dot local term modified by the presence of 
the magnetic field and the electrostatic interaction reads: 

Hd = ^ Sda^'ida- + Uoondtndi + Hpoi . (34) 

Here Sda = ed + (^e + ct^bB, (cr = + for f, f = — for 
4,) and is the Bohr magneton. Note that we included 
the Zeeman splitting in the spectrum of the dot, but we 
don't need to include it in the lead spectrum, since the 
density of states of metallic leads can be considered as 
constant at the Fermi energy. The last term describes 
tunneling and generalizes the one entering Eq. (|22p to 
take into account the spin: 

Ht^ ^{<.exp[*^(6„ + 6t)0,„]4^d, + h.c.}. (35) 

In the absence of magneto-elastic coupling, the Hamil- 
tonian ([5^ was introduced in Refs. to describe 

transport through a quantum dot in presence of a lo- 
cal electron-phonon coupling and later studied by several 
authors^i^>21^>23 



C. Combined coupling constant 



where a indicate the spin projection opposite to a. In 
Eq. (|38p we have introduced the complex electron-phonon 
coupling constants: 

aun = K + i^i/n ■ (39) 

Let us first discuss the effect of the prefactor 
exp(— 2i 0i,„A„ntiCT)- Since we are considering the 
limit where double occupied states are not accessible, it 
is easy to show that this term gives always 1 when evalu- 
ated with da on the three dot's basis states {|0), | f), | i)}- 
We will thus drop it in the following. The phase factor 
exp(— i0,y„A„) can be included in the definition of 
the tunneling amplitude t^. 

Combining the above results, we finally obtain the 
Hamiltonian 

H = ^OJnbibn+ ^ £,,ykcl^aC^ka +^£danda 
n v,k,(7 u 

+ Uoond\nd^ + Ht ■ (40) 
with the tunneling Hamiltonian: 

Ht^YI cxp[^(a^„6„ - alJ,l)\clf^^da + h.c.}. 

u.k.cT n 

(41) 

In Eq. (|1T|) the polaronic and magneto-elastic cou- 
plings enter on the same footing through constants a^m 
see Eq. ([39l). 



The expression for the Hamiltonian can be further ma- 
nipulated in such a way that the dimensionless coupling 
constants t/)„ and A„ introduced in subsections III Al and 
IIIBl respectively, can be treated on the same footing. 

For this purpose it is convenient to eliminate the linear 
polaron coupling in the Hamiltonian p4p by using the 
Lang-Firsov unitary transformation H UHU^: 

The transformation shifts the bosonic fields {UhnU^ ~ 
bn - Knd) giving for i/string + Hd'. 

Hd + -^string ^ ^nbliK + X] ^danda + Uoondtndi , 

(37) 

with £dcr = Sda - J2n ^n^n and Uoo=Uoo-'^ ^n^n- 

Both the energy of the dot level and the effective interac- 
tion between two electrons on the dot are renormalized. 
Note that, for Uoo < '^J^n'^n^^^ ^^'^ effective interac- 
tion can become attractive f^^^ Focusing our analysis to 
strong Coulomb repulsion regime we are not concerned 
by this possibility. 

Using Ud^W = d^e^" ^-('•'•-''n) and Baker-Hausdorff 
formula, the tunneling term in Eq. psp becomes 

^p-2.i:„ A„0.„„,, cxp[X(a,„6„ - Knbi)]da , (38) 

n 



D. Discussion 

Remarkably, the ratio between the magneto-elastic and 
polaronic coupling constants, cf. Eqs. and ([HT]) . is 
independent of the length of the nanotube. A crude es- 
timate shows that, for the fundamental bending mode, 
(j)o/^o ^ \/ni/m*{e/Cg\Vg\)BhaB/^o, where gb is the 
Bohr radius and m* is the nucleon mass. (To get it we 
used p ~ m* /ub, I ~ R± ^ and E ^ Ry/aj^, where 
Ry is the Rydberg energy.) In the ratio (j)o/^o, the small- 
ness of m/m* is compensated by the largeness of h/aB- 
Thus, magneto-elastic and polaronic coupling constants 
of the same order of magnitude could be realized, while 
being tunable independently through the magnetic field 
and gate voltage. It is thus interesting to investigate their 
interplay in the transport properties of the device. 

Let us now estimate the coupling constants quantita- 
tively. Assuming a linear mass density p of the order of 
10^^^ kg/m, corresponding to six carbon atoms/A, and a 
length L of 1 ^m, we find that the typical mass of a single- 
wall carbon nanotube is 10~^^ kg. The zero point motion 
Xno for a mode of frequency ujn/'^'^ = 500 MHz is 4 pm. 
For B = 40r, Eq. ^ yields 0o « 0.1 for the funda- 
mental bending mode [for which J^^ (ds / L) fo{s) ~ 0.83], 
while, for h = 200nm and CgVg — 5e, one finds from 
Eq. ([5T|) that Ag ~ 0.1. The magnitude of the magneto- 
elastic and polaronic effects is thus the same for these 
parameters. 
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It is probably useful to point out that the current 
dependence generated by the strong magnetic field dis- 
cussed in this paper and in previous works^^^— is totally 
unrelated to the Aharonov-Bohm effect discussed in the 
literature in several papers (see Ref. [s^ for a review). 
The effect discussed there is due to the interference of 
the electrons in the Aharonov-Bohm loops generated by 
the carbon crystallographic structure. The presence of 
phonon modes is not required for this effect to be present, 
while it is crucial for the magnetoresistance discussed 
here. It should be also noted that the two contributions 
are of course present at the same time in a realistic de- 
scription of the carbon nanotube, but the two effects can 
be distinguished experimentally since they have totally 
different thermal and voltage dependences.— 

Concerning the model for the description of the nan- 
otube deformation, in this paper wc assumed that the 
standard theory of clastic continuous media applies. This 
should be correct, at least for not too short nanotubes, 
where an atomic description becomes necessary. We only 
considered the case of a rigid rod with doubly clamped 
boundary conditions. It is well know that it is possible 
to observe a crossover from this regime to the regime of 
tense stringi^ The effect of a residual stress on the string 
can be taken into account by modifying the Hamiltonian 
(jl]).— The structure Eq. ([22|) remains valid in this case. 
The effect of the tension modifies only the profile func- 
tion /,i(s) entering Eq. for the eigenmodes. The 
elongation of the nanotube induced by the external gate 
and magnetic field, with typical amplitude ^ XnXno and 
4>nXno respectively (we took / ^ ea;„ for the last esti- 
mate), also generates an internal stress. The nonlinear 
effects induced by this stress are safely neglected when 
these amplitudes don't exceed i?_L,— which is the case 
for the parameters given above.. 

The amplitudes of the effective phases in the tunneling 
part of Eq. (j22p only coincide for symmetric boundary 
conditions. For instance, in the case of a nanotube that 
is clamped on one side and hanging on the other side, no 
specific relation would exist between the phases in the 
left and right terms. 

Surprisingly the Hamiltonian (|40 )) - (|4T|) with A„ = 
has been already introduced before (in the spinlcss case) 
for a quite different problem. Imam et al^ considered 
resonant tunneling in presence of an electromagnetic en- 
vironment associated with the fluctuations of the bias 
voltage, and described by a bath of phonon modes. Nev- 
ertheless, the effect of the environment on the conduc- 
tance was only addressed at large bias voltage or far from 
the resonance in that work. Moreover in this paper we 
address the interplay of the magnetic and polaronic cou- 
pling, that was not present there. 



(Note that tildes that appear there will be omitted from 
now on.) Here we present the general approach that will 
be used in the next Sections. 

The problem at hand is a generalization of the well 
studied polaron problem. It has not been solved exactly 
yet, to the best of our knowledge. Nevertheless, it is pos- 
sible to obtain approximate results in different regimes. 
When the coupling constant is small, a perturbativc ex- 
pansion is possible; it has been performed both for the 
polaronioi^ and the magneto-elastic cases.— In principle 
the method used in those previous publications can be 
extended to investigate the case when both couplings are 
present. However it works only when the predicted effects 
are small and we will not use it in this paper. Alterna- 
tively, the dynamics of a single electron in the presence 
of the phonon bath can be obtainedj^i^ However, the 
results obtained for the phonon-assisted resonant tunnel- 
ing can only be applied far from the resonance, when 
effects associated with the Fermi sea in the leads can be 
neglected. Another tractable limit is the case of high- 
temperature incoherent tunneling limit (T 3> F), which 
has been discussed in Ref. [13 for the polaron problem 
and does not require that the electron-phonon coupling 
is perturbative. Here, we generalize this approach to de- 
scribe the effect of the interplay of the two couplings. 

In this incoherent regime the electronic environment 
at temperature T has the time to suppress the coherence 
of the quantum evolution between tunneling events. The 
equations of motion for the off-diagonal elements of the 
density matrix in the basis of the eigenvectors of Hq — 
H — Ht decouple from those for the diagonal elements. 
It is then possible to obtain the equations of motion for 
the diagonal elements {i.e. the occupation probabilities) 
alone. 

The rates can be obtained by standard perturbation 
theory. The time evolution of any eigenstate |i) of Hq 
can be written in terms of the resolvent: 



m) 



dE 



-iEt 



2tt E - Hq - Ht + ill 



(42) 



The probability that the system is in another state |/) 
after a time t is 



P, 



dE 
2^' 



-iEt 



{f\G"{E) + G'{E)T{E)G'{E)\t) 



(43) 

where T{E) ^ Ht + G^{E)HtG°{E) + ... is the T- 
matrix and G^{E) ~ {E — Hq + i?/)~^ is the free retarded 
propagator. Neglecting the E dependence of T{E) in 
Eq. gives the standard Fermi golden rule with T{Ei) 
replacing Ht, that is, Pif = tWif with 



III. FORMALISM 

Our aim is to derive the current-voltage characteristics 
for the problem described with Hamiltonian (P(7|) - (PT]) . 



W,f=2n\T,jiE,)\^SiE.,-Ef)., 



(44) 



and Ei the eigenvalue of Hq related to the state i. By 
inserting the resolution of the identity in terms of the 
eigenstates of Hq, one can then write for the first two 
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orders: 



{HTU{HT)nf 



Ej — En 



(45) 



These two terms give the tunnchng and cotunnehng con- 
tributions to the transition rate, respectively. We wih 
investigate the effect of the tunnchng term in Section IIVI 
and then the corrections given by the cotunnehng terms 
in Section [V] A comment is at order. The expression 
dU]) holds if the frequency dependence of T(a;) is suffi- 
ciently smooth. This is not true if one intermediate state 
has the same energy as the final state. The problem has 
been discussed in Ref. [23 and it requires a proper renor- 
malization of the divergences of the integrals appearing 
in the calculation of the transition rates. The regulariza- 
tion for the transition rates associated with cotunnehng 
processes is done in Appendix IB] 

The system state is thus completely described by the 
probability of being in one of the eigenstate of Hq. We 
will assume that the lead remains in thermal equilibrium 
during the evolution, thus the only relevant degrees of 
freedom left are the dot and oscillator degree of free- 
dom. In the following we will focus on a single phononic 
mode, for which we define lo the resonance frequency, (p 
the magneto-elastic coupling, and A the polaronic cou- 
pling. In order to fully describe the system we finally 
define the probability Psn (t) that the dot is in one of the 
three available states (|s) = |0) , jj,) , |t)) and the oscilla- 
tor in the phonon state n (|n) — b^"' \0) / Vn\) . The rate 
equations read: 



dt 



- i I Psn - £ PsrrA ■ (46) 

^ \ m=0 J 

Here Wsn,s'm is the transition rate from the electronic 
dot state s and phonon state n to another configuration 
(s', m). Following Ref.[27l a phenomenological relaxation 
term has been introduced in Eq. to describe the 
coupling of the oscillator with the non-electronic envi- 
ronment at the origin of an intrinsic finite quality factor 
= wT/27r for the oscillator. The function 

Pn gives 

the equilibrium phonon distribution function: 

Pj^'i^ = (1 - e-"/^)e-""/'^ . (47) 

Given the rates, equations ([^S]) can be solved; the aver- 
age current is then obtained from the stationary solution 
for the probabilities. In the next two sections, we discuss 
the result for the current at the sequential (Sec. lIVp and 
cotunnehng (Sec. |V| orders. 

IV. SEQUENTIAL TUNNELING REGIME 

The sequential tunneling regime corresponds to keep- 
ing only the lowest order term in the expansion (|45p that 



determines the transition rates. In this regime we find 
that the current-voltage characteristics depends on a sin- 
gle coupling constant which does not allow to distinguish 
the polaronic and magneto-elastic effect from each other. 



A. General formula 

The explicit form of the transition rates in the sequen- 
tial tunneling limit reads: 



sn.s' m 



sn.s' m 



(48) 



with 



x6{Ef^ - E,_^ + Eds' - £ds + uj[m - n]) . 

(49) 

Here v and f^, are initial and final fermionic states in 
lead V, with energies i?;,^ and Ef^^, P^'^'^^i^) is the equi- 
librium distribution of state i^,, and for convenience we 
introduced e^o = 0. The non vanishing rates take the 
form: 



W, 



On,crm 



riu {Sda + ^[m - n]) 



(50) 



am.On 



M'^^^ [1 - {sda + w[m - n])] , 



(51) 

where n^{e) = \j\e^'^~^^^l'^ + 1]"-'^ is the Fermi distribu- 



^aiyb—a*b' 



{n\e""" 'I'm), and we recall 
that = 2TTm'^N^. 

The absolute value of the matrix element Mnm depends 
only on |q!i^| = a, where 



(52) 



is the r.m.s. of the polaronic and magneto-elastic cou- 
pling constants. This can be readily be shown by absorb- 
ing the phase of in the operators b and &^ appearing in 
the definition of Mnm- This is a surprising results, which 
holds only at this sequential order in the tunneling expan- 
sion. For instance, it was shown in Ref. [2^ that there 
are qualitative differences between the magneto-elastic 
and the polaronic problems: at resonance, the maximum 
current depends on the asymmetry of the electronic cou- 
pling between the dot and the leads in the first case and 
does not in the second case. It is thus clear that to higher 
order differences have to appear. Technically this means 
that the current should depend also on the phase of the 
coupling ai/, as we will find at the cotunnehng level in 
Sec. El 

In the reminder of this Section we will study the evo- 
lution of the current- voltage characteristics as a function 



of the magnetic field. For simplicity we assume that the 
relaxation time introduced in Eq. (j46p is shorter than the 
typical tunneling rate, so that the distribution function 
can be written: 



Ps 



p p(eq) 



(53) 



The rate equation (|46p for the stationary solution 
{dpsn/dt = 0) can then be cast in the form WP = 0, 
where P = (Poi-Pfi^i) (with the normalization condi- 
tionPo + n + n = l)> 



/ 



V 



to 



w, 



to 



to 




and we define the reduced transition rates: 

ss' / J n "sn,s'm^ 



(54) 
(55) 



By performing the sum over the phonon states the rate 
take a form similar to the non interacting cascj^^ ^i)^ — 
fi,{eda) and W^^^' = w^^J e^^^'-^^^^^ with 



n^(e) ^ j dE V{E) n^{e + E) . 
s 

dte}^\A,(t)Al{Q)) , 



Here the function V reads 
1 



V{E) 



27r 



(56) 



(57) 



with A^{t) = (f^uh[i)-aih\t) _ Fourier space, V{E) 
takes the form of a a sum of delta functions centered 
at integers multiples of the oscillator frequency oj. In 
the limit T ^ w, on which we will focus hereafter, the 
oscillator is in the ground state and 



/T/ 

V{E) ^ E e-" — -6{E - mw) 



(58) 



The current, say for instance at the left lead, can be 
written in terms of the probabilities and the tunneling 
rates: 

oo 

^ = e XI H (^On.amPOn " W^i2j,0nP<Tm) ■ (59) 
(T=t,t 

In the case of interest here of the equilibrated oscillator 
one obtains: 



1= E k^-p^)i^ 

(T=t,t 



(60) 



with the current for each spin-channel: 



(61) 



and the stationary probabilities: 



P. 



(ed5-;V')/T 



T,^.,.' r^r^/ni.(edt)n^'(edt) [nu\£dt)n^,^{£di) - l] 

(62) 

with Pq = 1 ^ -Pt ^ -^t- Formula (|60| for the sequential 
current has a simple interpretation: the current associ- 
ated to the tunneling through one spin state and in ab- 
sence of any correlation between the electrons is reduced 
by a factor which is just the probability of occupancy of 
the dot with opposite spin. Indeed, Icr corresponds ex- 
actly to the sequential current for the spinless problem 
in presence of electron-phonon interaction. 

From the previous expressions, we can obtain the linear 
conductance in the sequential tunneling regime: 



G — Gq 



T 1 + 2e-(-<i-'')/'^ coah{iiBB/T) 



(63) 



where Go = e^FiF^/r^, Sd = (e^t + edt)/2 and n = 
(/i; -I- iir)/'^- Eqs. are the central result of this 

Section. In the following we discuss the current in differ- 
ent regimes. 



B. Linear conductance 

The linear conductance generically displays a reso- 
nance as the energy of the discrete level in the dot 
is varied, e.g. with a gate voltage. 

In the absence of electron-phonon interaction;^ the 
conductance is peaked around £d~ n — {T/2) log 2, with 
width ^ T and maximum height Gmax — 2Go(F/T)/(3 + 
2V2) at zero magnetic field. It is peaked around around 
Sd — fJ- — (J-bB (assuming _B > 0), with width ~ T and 
maximum height G^ax = Gor/(4T) at large magnetic 
field /^B^ > r. 

The main effect of electron-phonon interactions is to 
suppress the peak height: At low temperature (T ^ w), 
we find with help of Eqs. ([55]) . ([55]) that the suppression 
factor is e~" 



C. Non-linear IV characteristics 

In this section, we assume that the bias voltage is ap- 
plied symmetrically, with fii^r = A* i eV/2. 

In the absence of electron-phonon interaction and at 
zero magnetic field, resonance lines at £d — M — 
with width ~ T, separate the regions where the current 
is blocked and those where it reaches a finite value / = 
2eF,Fj(F-FFz) or / = -2eFiF^/(F -h F^), depending on 
the sign of the bias voltage^i At finite magnetic field, 
the regions where the current is blocked and those where 
down spins only can flow through the dot are separated 
by the lines e^— /i— /^gi? ~ ±eV/2. Inside the conducting 
regions it appear bands where an additional conduction 
channel opens for up spins. At 1^ > (respectively V < 
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0), these bands stand between the hnes — n~ eV/2 = 
±fiBB {sd — H + eV/2 = ±^bB); the current reaches 
/ er,r,./r (/ = -eViTr/r) inside the bands. The 
contrast between the different lines in tlie gate and bias 
voltage dependence of the differential conductance, see 
Fig. [21 is readily explained by the amplitudes for the 
current found above between the lines. 



( a//av )/Go 

0.04 




2 4 6 

(8d-^)/T 



8 10 



FIG. 2: Differential conductance in the absence of electron- 
phonon coupling {a = 0) as a function of gate and bias volt- 
ages in the sequential tunneling regime, for fisB/T — 5, 
Fi/T = F,./T = 0.1. 



The case of finite electron-phonon coupling and zero 
magnetic field was discussed in details in RefsHHH. 
Note that Equation ([6T|) for the current simplifies into 



2eF,F. 



ni{£d)ni{ed) (e 



{ed-tJ.R)/kBT _ piea-fJ.L)/kBT\ 



(64) 

The main features of the current are (i) the presence of 
inelastic lines (vibrational sidebands) at Ed — fJ- + nuj = 
±eV^/2 in the conducting regions and (ii) a suppression of 
the current at the low bias voltage (Frank-Condon block- 
ade) when the electron-phonon coupling is large enough. 
Fig. 13] illustrates point (i) for intermediate electron- 
phonon coupling. 

We can now understand the general features of the cur- 
rent when both the electron-phonon coupling and mag- 
netic field are present. They look different depending on 
the ratio fisB/ui. When uj > fisB, the vibrational side- 
bands are split due to the removal of the spin degeneracy 
of the single level, see Fig. 01 When fisB > lu, vibra- 
tional sidebands are now present as replicas of the two 
main resonance lines, see Fig. |5l 

In conclusion of this Section, we found that the po- 
laronic and magneto-elastic couplings cannot be distin- 
guished in the sequential tunneling regime, apart from 
the removal of the spin degeneracy induced by the Zee- 
man effect. In the next Section, we will obtain a different 
conclusion when higher order processes in the tunneling 
expansion are considered. 
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FIG. 3: Differential conductance at zero magnetic field as a 
function of gate and bias voltages in the sequential tunneling 
regime, for a = 1, ui/T = 20, and Fi/T = Fr/T = 0.1. 
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FIG. 4: Differential conductance as a function of gate and 
bias voltages in the sequential tunneling regime, for A = 1, 
= 0.4, hbB/lo = 0.2, Lo/T = 40, and Vi/T = Tr/T = 0.1. 



V. COTUNNELING 

We first discuss the general approach to consider the 
cotunneling corrections to the rate equations. Then we 
show that, in a specific regime, the cotunneling process 
dominates the total current. Finally, we concentrate on 
the elastic tunnel current to clarify the differences be- 
tween the magneto-elastic and polaronic interactions. 



A. General formula 

The explicit form of the transition rates induced by 
cotunneling processes is given in Appendix |B1 

Assuming that the phonons are equilibrated, we find 
that the rate equation (|46p including cotunneling events 
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I 



due to the electron transfer via cotunneling: 

+ e E (69) 

and the effect of the variation of the occupation proba- 
bilities produced by the cotunneling processes: 

<5/2 - e E {W^'^SP, - (70) 

cr 

In the next Section, we discuss the results for the current. 



FIG. 5: DifTerential conductance as a function of gate and 
bias voltages in the sequential tunneling regime, for A = 1, 
(f) = 0.4, fj.BB/uj = 2.75, uj/T = 20, and Ti/T = T^/T = 0.1. 



reads 



(w + m)(p + 5P) = 0, 



(65) 



where SP is a correction to the occupation probabilities 
P found in Sec. IV] (with SPo + 5P^ + SP^ = 0) and 



SW 












wi^^"^ -wl^""^ 



(66) 



and we introduced the reduced cotunneling rates 



sn.s'm ' 



(67) 



corresponding to the transmission of one electron from 
lead v to lead v' without variation of the dot state, s = 
s' = 0,t,i, or with a spin-flip process between a ='f and 
a' =i 

Considering that the cotunneling rates are smaller than 
the sequential tunneling rates by a factor ^ T/T^ we ne- 
glect a term (5W(5P in Eq. (|65l) and we find that the 
correction to the occupation probabilities solves W(5P ~ 
— (5WP. Note that, when the oscillator is at thermal 
equilibrium, the cotunneling processes affect the steady 
occupation probabilities P only via the spin-flip pro- 
cesses. In the limit of vanishing magnetic field, we have. 



by symmetry, W^.||"^'' 



and Pf = Pi . Then the co- 



tunneling rates completely cancel from Eq. (|65p , in agree- 
ment with the results of Ref. [l^. 

According to the previous discussion, the cotunneling 
correction to the current can be decomposed in two parts: 

SI = 5Ii + Sh ■ (68) 

They represent respectively the correction to the current 



B. Discussion 

The above formalism allows us to obtain the current 
including both sequential tunneling and cotunneling con- 
tributions, for arbitrary parameters of the junction. An 
example is shown in Fig. [51 In general, cotunneling 
processes produce a smoothing of the features already 
discussed in the sequential current, that is the steps as- 
sociated with the vibrational sidebands and the Zeeman 
splitting. However they dominate the current in the re- 
gions where it is strongly suppressed at the sequential 
level, see inset of Fig. [51 

In the following, we will concentrate on the region 
where the cotunneling terms dominate the current. To 
simplify the expressions for the current, we will further 
assume that the gate voltage is tuned in the region where 
the dot is mostly empty (Pq « 0) and the magnetic field 
is sufficiently large that only the down spin state for 
the singly occupied dot can be accessed virtually, that 
is when — ji, fisB 3> P, eV. In addition, we will con- 
sider that T ^ w so that the oscillator is in its ground 
state. Then the expression for the current (noted / 




FIG. 6: Bias voltage dependence of the current with (straight 
line) and without (dashed line) the cotunneling contribution, 
for A = 1, = 0.4, Ti/T = r,./T = 0.4, lu/T = 15, {sa - 
fi)/ui = 5, and hbB/u) — 2. Inset: Zoom of the low- voltage 
region in log-scale. 
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instead oi SI from now) simplifies to 

/ = e (M/^r - , (71) 

with 

^(i/i/') _ i-ui^i^' sr^ M^"'^* M^''"^ M^"'^ m'^''"^* 

^^00 — ^ -'"mg -""Op -'"mp 

m,p,g— 

xT{n^;H„' +mui;Sdi + quj;Sdi+pi^)- 

(72) 

where the function J-" is defined in Appendix B. A finite 
electron-phonon interaction results in the suppression of 
the cotunneling current, compared to the current in the 
absence of that interaction 7^") = (/or/47r2r){*'[l/2 + 
«(/ii - £di)/(27rT)]-*'[l/2+i(/x, - edi)/{2TTT)])}, where 
Iq = eTiTr/T. Similarly to the results found in Sec |IV] 
we find that the polaronic and magneto-elastic couplings 
cannot be distinguished close to the charge degeneracy 
point, at — jj, <^ uj. Focussing for instance on the 
elastic contribution to the cotunneling current, 

/ = /^lA^^plMMp = /(0)e-2a^, (73) 

we find that the suppression factor is the square of the 
one found in the sequential tunneling regime. This ef- 
fect can be understood by (i) projecting Hamiltonian 
(P(I l) - (|4ip over the oscillator's ground state^^ thus ob- 
taining an effective Hamiltonian for a non-interacting 
resonant level with renormalized tunneling amplitudes 
t'^ = exp (— a2/2), and (ii) recalling that the current 
is obtained as a fourth-order efi'ect in these tunneling am- 
plitudes. 

Instead, far away from the charge resonance at e^j, — 
/i ^ cj, we find 

00 

q,p=Q 

= /(•')|(0|A,4|0)|2 = /(o)e-20' (74) 

(here A^, = g'^fb-oi^b'' -j^ -ftrhich is only suppressed by the 
magneto-elastic coupling. The absence of an exponential 
suppression of the cotunneling current in the polaronic 
case was noted in Ref. [l^- The magneto-elastic effect can 
be understood by accounting for the virtual transitions 
to the singly occupied states of the dot perturbativelyi^ 
In general, the polaronic and magneto-elastic couplings 
suppress the cotunneling current in different ways at not 
too large oscillator's frequency. The difference is traced 
back to the way the couplings appear in the tunnel- 
ing Hamiltonian (|^T|) . In the polaronic case, an elec- 
tron crossing the junction is subject to opposite phases 
±iA„(6„ — 6jj) at the tunnel barriers with the leads and 
the dot, which can result in a compensation of their effect. 
Instead, in the magneto-elastic case the electron feels the 
same phase (l){b + b^) at each barrier. Indeed, for the 
same amplitude of the coupling constants, we find that 



the suppression of the cotunneling current is stronger for 
the magneto-elastic interaction than for the polaronic in- 
teraction, see Fig. [T] 




eV/co 



FIG. 7: Bias voltage dependence of the current for different 
electron-phonon couplings and Fj/T = Tr/T — 0.4, oj/T — 
15, {sd — fJ-)/^ = 5, and ^bB /lo — 2. 



VI. CONCLUSIONS 

In conclusion, we have studied how the quantum vibra- 
tions of a suspended carbon nanotube in the quantum dot 
regime affect the charge transport through it. We have 
paid special attention to the respective signatures of the 
polaronic effect that arises from the position-dependent 
capacitive coupling of the nanotube with a nearby gate 
and the magneto-elastic coupling in the presence of a 
transverse magnetic field. While both couplings act sim- 
ilarly in the suppression of the current at the sequential 
tunneling order, their effect is qualitatively different in 
the cotunneling regime. 

Our estimates for the polaron and magneto-elastic cou- 
pling constants show that these effects may be relevant 
to determine the current-voltage characteristics through 
suspended nanotubes. Thus, charge transport could be 
efficiently used to demonstrate the quantum nature of 
the mechanical vibrations of these tiny objects. 
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Appendix A: Rigorous derivation of Eq. ([S]) 

In this Appendix, the effectively one-dimensional 
Hamiltonian ([S]) for electrons in a suspended nanowire 
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is derived and its conditions of validity are obtained. 

For this, we start again with Eq. ([1]) and then per- 
form a unitary transformation T-L — > e^^7ie~^^ where 
S — Pyu{x) — eBU{x), to obtain 

^ ^ [Px-VydMx)-eByY+p\ ^ ^^^^^ 



2m 



E 



[Pn-Pyfn{x)+eBFn{x)Y , MuolXl 



2M 



(Al) 



Here p± — {py,pz). Assuming a strong confinement in 
the nanowire, we project Eq. (|A1I) on the ground state 
for the electron's transverse motion, and we obtain 



n 



[Px - eB{y)o] 



2m 

[P„ + eBF„(x)]2 



2M 



+ 



+ 



e5({2/>Py})o 



2m 



{pI)o 
2M 



{pI)o 
2m 



with 



£t — £o 



e'B' 
2m 



[{y')o {y)l] 



(A2) 



(A3) 



In Eqs. (|A2|) - (jA3|) . we denoted (a)o = / dr ±XoO-Xo for an 
arbitrary function a(r±,p±), where r± = {y,z) and xo 
is the ground-state wavefunction for transverse motion 
that solves the Schrodinger equation: 



+ t/conf (r±) 

Zm 



(A4) 



(Note that (py)o ==0.) The first term in Eq. ([M]) is the 
ground-state energy Eq for transverse motion, the second 
term is a global diamagnetic shift. Both will be absorbed 
into the chemical potential from now. 



We now argue that the terms in the two last line of 
Eq. (|Aip can be safely neglected: 

The first term would lead to a shift in the equilibrium 
position of the nanotube that is characterized by a rela- 
tive displacement SXn/Xno ~ {M/m){Xno/L)'^4'n of the 
eigenmode n. We also anticipated that the relevant mag- 
netic field dependance in the model enters through flux 
~ BLXno, see Eq. 

The second term (like the third one) acts even in 
the absence of magnetic field. It generates an energy 
shift Ae when the electron is in the suspended sec- 
tion of the nanowire, with typical amplitude Ae/cj„ ~ 
{Xno/L±)^{L/a). Here, L± is the transverse width of 
the wire and a is a cut-off at short distances (on the 
nanoscale) which regularizes the sum J2n fni^) ~ ^(^) 
that diverges formally in the approximation of continu- 
ous medium leading to Eq. The parameter a repre- 
sents the microscopic atomic distance, the natural length 
scale at which the continuous model ceases to be valid. 

The last term leads to a renormalisation of the oscilla- 
tor's frequency when the electron is in the suspended 
nanowire, with typical relative magnitude dujn/uJn ~ 
(Af/m).(X„o/i)^.(^„o/^±)^- We can check that all 
terms give small correction in the case of the single- wall 
carbon nanotube with parameters given in Sec. IIIDI 

Taking into account these simplifications in Eq. (|A1[) 
and making another unitary transformation Ji — > 
e^^ 'He~^^' , with S' = eBx{y)o, we finally recover Equa- 
tion (E). 



Appendix B: transition rates associated with 
cotunneling events 

In this Appendix we obtain the transition rates as- 
sociated with cotunneling events following the method 
discussed in Ref. [13 

In the absence of sequential tunneling term in Eq. (|45l) , 
the transition rate for the dot-phonon state to evolve 
from (s, n) to (s', m) while an electron is transferred from 
lead ly to lead i^' is 



wiZL^"^^ E E P^"''Hiu)P^"''H^^')\{UJ^':s\n\HTGoHT\i.,i.',s,m}\^S{AE,+AE,,+Ass'+A,nn) 

(Bl) 

where AE„ = Ef^ - Et^ , A^n = w[m - n], and Ass> Eds - £ds'- 
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The non- vanishing rates take the form 



w. 



On,Om 



an.a' 7n 



^ j den.(C)[l-n,,(e + A„„)] 



E 



27r 



den,(C)[l-n,,(e + A 



E 

9=0 



A 



■mq 



11] 



(B2) 



(B3) 



These expressions diverge as 77 — >■ 0. A regularization scheme to remove this divergence while the sequential tunneling 
events are properly taken was discussed in Ref. [l^- It consists in removing an C'(l/?7)-term, which corresponds to 
sequential tunneling events and is already accounted for in the rates (j49|) . while keeping the next order term [0(1) 
term]. On the end, the regularized cotunncling transition rates read 



where 



W, 



r.r 



On,Om 



2^ / , / , "'^nq ""-mq ^^^np -^"mp iPi^jpi^' '-^nnn^da '-^nqi^da '-^npj 



r r / 

27r 



QO 



EM^'''>*M^'^'^M^''^M^'''^*F(u ■ ;y / - A /-A ■ / + A ■ , + A ) 



p,q=0 



T{Ei, £2,61,62) = limRe 

r|—^0 



dE 



i{E - Ei)[l - n{E - E2)] 



{E-ei+irj){E 



if]) 



-oi- 

V 



nB{E2-E^) 



Re 



£1 



, , 1 E2-ei 

* 

2 2ttT 



, , 1 E2-e2 

- * 

2 2ttT 



, , 1 Ei-ei 

- + I— 

2 2-kT 



2 27rT 



UsiE) is the Bose distribution function, and ^ is the digamma function. 
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